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Abstract 



We consider the production of muon pairs from the scattering of pious on 



Oh! 
> 

X 

■ longitudinahy polarized protons. We calculate the cross section and the sin- 

gle spin asymmetry for this process, taking into account pion bound state 
effects. We work in the kinematic region where the photon has a large longi- 
tudinal momentum fraction, which allows us to treat the bound state problem 
perturbatively. Our predictions are directly proportional to the pion distribu- 
tion amplitude. A measurement of the polarized Drell-Yan cross section thus 
allows the determination of the shape of the pion distribution amplitude. 
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I. INTRODUCTION 



Spin effects are known to provide sensitive tests for the underlying particle theories in 
many cases. One may recall the spin asymmetries in pion-nucleon scattering which cannot 
be described in the simplest pole approximation of the Regge theory and require to take into 
account cuts. In the case of QCD the so-called EMC Spin Crisis was a signal of a breakdown 
of the naive parton model and of the importance of the incorporation of such a subtle field 
theory effect as the axial anomaly 

A less popular, but also extremely interesting example is the problem of single transverse 
spin asymmetries in reactions mediated by strong interactions. As these effects require a 
mass parameter and an imaginary part in the scattering amplitude, they are proportional 
to rrigas in the 'naive' perturbative approach. However, a careful application of QCD fac- 
torization at twist-3 level removes both of these small parameters 0. The collective gluon 
field of the hadron is shifting the current quark mass nig to a mass parameter of the order of 
the hadron mass and, simultaneously, provides the imaginary part. Technically, the latter 
appears when the pole of the quark (or gluon 0) propagator is integrated over the light-cone 
momentum fraction. 

A similar imaginary part is present in the QCD higher twist contributions to dilepton 
production in pion-nucleon scattering, due to the integration over the light-cone momentum 
fraction of the quark in the pion. This effect leads — in the case of unpolarized nucleon 
targets — to a significant contribution to the dilepton angular distribution [§], which is 
sensitive to different ansatze for the pion distribution function. 

It is interesting to look for single spin asymmetries for which the imaginary part does not 
constitute a (large) correction, but the whole contribution. The natural candidate is the spin 
asymmetry in the production of dileptons from the scattering of pions on a longitudinally 
polarized nucleon. This spin asymmetry has been thoroughly studied in perturbative QCD 
I^J^. The analogous mesonic processes have been discussed in 0. In the present work we 
study this single spin asymmetry, taking into account pion bound state effects. 
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II. CALCULATION OF THE HADRONIC TENSOR 



We will consider the inclusive production of muon pairs from the scattering of a pion on 
a polarized proton (the Drell-Yan process), 

7r-(P^) +p(Pp, se) ^ 7*(Q, e)+X-. + /^-(g.) + X, (2.1) 

where Si denotes the degree of longitudinal polarization of the proton, and all momenta refer 
to the overall cm. system. It is convenient to write the differential cross section for ( |2.1| ) in 
terms of a hadronic and a leptonic tensor, 

where s is the hadronic cm. energy. The leptonic tensor has the well-known form 



If^u ^ j^^/.^. ^ _ ^m-q2/2) . (2.3) 

Note that L^'^ is symmetric, since we only consider the exchange of a virtual photon (not a 
Z boson) and since we sum over the polarizations of /i^. 
The hadronic tensor W'^'^ may be written as 

X 

{^{Q,e^),X{Px)\T\n-{P^),p{P„se)) 

{^iQ,en,XiPx)\T\n-iP^),piP„s,)y. (2.4) 



Only the symmetric part of W'^'^ is relevant after contraction with L^^, of eq. ( p.3|) . Therefore 
W'^'^ in the following is implicitly assumed to be symmetrized. 
The angular distribution of in (|2.1| ) may be parametrized as 

^'^ oc 1 + A cos^ 6' + /i sin 29 cos H — sin^ 9 cos 20 



dQ'^dQ'^dx Ld cos 

sin2^sin0 + — sin^ 6'sin20. (2-5) 

Here, 9 and are angles defined in the muon pair rest frame, and A, /i, u, fi, v are 
angle independent coefficients. They depend on (the virtuality of the photon), Q\ (the 



squared transverse momentum of the photon in the hadronic c.m.s.), and on xl = 2Ql/^/s 
(the longitudinal momentum fraction of the photon). We take > for a photon moving 
forward with respect to the pion in the hadronic cm. system. The normalization of the 
cross section can be determined from 

Q^da _ (3 + A)iV 

dQMQldxL 96(2vr)4 + Ql + xls/A ^ ' ' 

The normalization N and the angular coefficients are related to the helicity amplitudes 
of the hadronic tensor: 



N = 


Wt + Wl = 


-- e^(+l)PF^^e:(+l) + e^(0)iy'^%(0) 


X = 


N- 


-1 {Wt - 


Wl) = iV-i{e^(+l)iy^'^e:(+l) - e^(0)iy^%(0)} 


= 


N- 


''Wlt = 


(v^iV)-i{e^(0)iy'^%(+l) + e^(0)iy^'^e:(+l)} 


u = 


N~ 


^^Wtt = 


N-'{e,{+l)W^Xi-^) + e;(+l)W^^%(-l)} 


= 


N- 


-'Wlt = 


z(v^iV)-i{e^(0)iy'^%(+l) - e^(0)iy^^e:(+l)} 


v = 




'^Wtt = 


^N~'{e,{+l)W^''e:i-l) - e;(+l)H^'^%(-l)}, 



(2.7) 

where the polarization vectors e'^(±l) = (0, =Fex + iey) and e^(0) = (0,6^) are determined 
by specifying the coordinate axes ej {i = x, y, z) in the muon rest frame. 

The coefficients /2 and v are nonzero only in the polarized Drell-Yan process and are 
induced by absorptive parts in the scattering amplitude T. In the parton model, the leading 
contributions to /i and v come from the interference of Born diagrams with the absorptive 
parts of one- loop diagrams 0. The relation of the single spin asymmetry 

^ _ d(y{si = +1) - dajse = -1) ,^ g, 

~ da{se = +1) + da{se = -1) ' ^ 

to the angular coefficients is simply given by 

^ /i(s£ = +1) sin26'sin0 + iz/(s£ = +1) sin^ 6'sin20 (2 9) 

1 + A cos^ + fi sin 26 cos + |z/ sin^ 9 cos 20 

For the unpolarized Drell-Yan process the parton model predictions |^ fail to describe 
the data on A, fi and 0. In contrast, taking into account pion bound state effects 
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yields reasonable fits to the data An interesting feature of this bound state model is a 
nonvanishing absorptive part in leading order which is proportional to the pion distribution 
amplitude ip. This property directly leads to nonvanishing ft and u which are strongly 
sensitive to if. In the following we will shortly describe the salient features of the bound 



state model for the pion. Further details may be found, e.g., in |T^, [Q. 

The bound state problem may be treated perturbatively if the momentum fraction of the 
quark from the pion is large, which is the case for large xl- The dominant contribution to 
the Drell-Yan process comes from the annihilation of a quark with an antiquark |Tl| . In the 



region of large xl, the diagrams of fig. la,b. then give the leading contributions. We will 
explain the formalism for the annihilation of a m quark from the proton with the u quark 
from the pion; the other possible annihilation process of a d quark from the proton with 
the d quark from the pion is treated in complete analogy and its contribution will later be 
added incoherently. 

In diagram la we see that the u quark propagator is far off-shell, p| = —Q^/{1 — Xu), 
if the light cone momentum fraction of the quark from the pion x^ ~ Xl is close to 1 | |10[| . 
Thus, the pion can be resolved by a single hard gluon exchange |]T^. The second diagram, 
fig. lb, is required by gauge invariance. 

The hadronic tensor W^" is obtained by a convolution of the partonic tensor w'^'^ with 
the corresponding parton distribution functions for the polarized protonQ. The partonic 
tensor w'^'^ is computed from the product M^M*^ . is the amplitude for the reaction 

M^ + TT- ^7*(g,e^)+X, (2.10) 

where the arrow indicates the polarization of the parton. This amplitude is obtained 
by convoluting the partonic amplitude T^{u^ + urf — > 7* + with the pion distribution 
amplitude ipi^z^Q"^) |T^, 



^Here we take into account the correct partonic flux factor which amounts to multiplying with 

S/S, S = {P^ +Puf- 

5 



^' = J^ ^ ) / ^ ) = 1' (2.11) 

4iVc JO JO 

where ~ Q\/{1 — Xu) is the cutoff for the integration over soft momenta in the definition 
of if, ~ 133 MeV is the pion decay constant, and A^^c = 3 is the number of colors. 



A detailed discussion of the scale setting problem for exclusive processes is given in |13 
Using the commensurate scale relations in the calculation of ( |2.11| ) would allow to eliminate 
ambiguities caused by the factorization and renormalization. However, this does not change 
the qualitative features of our predictions and thus it will not be discussed in this paper. 

The imaginary part in this model arises because the internal quark line of fig. lb can go 
on-shell. After integrating the hard scattering amplitude according to equation ( p. Ill) , 



we are left with a regular amplitude which has an imaginary part proportional to the 
pion distribution amplitude at the point where the quark propagator gets singular. 

We are now ready to present our analytic results. We relegate the full result for the 
hadronic tensor to the appendix, and confine ourselves here to a discussion of the angular 
coefficients /i and v. To present our results for the angular distribution we choose the 
Gottfried- Jackson frame where the axis is taken to be the pion direction in the muon rest 
frame and the e.x axis lies in the tt~P plane such that the proton momentum has a negative 
X component. From eqs. ( |2.7| ), and ( |A2| ) we get the angular coefficients 

—2se,pxFT[ip{x) 



[1 - xy [{F + Re I{x)y + 7rV(5)'] + p^x^F^{A + p^) 

X 



V = 2pp. (2.12) 
Here, p = Qt/Q and a; is a function of xl, Q'^/s and p 



xn lXL + Jxl + 4Q^s-\l + p^^ 



Furthermore, 



z 



m= Tdz , ^("'^'^ , (2.14) 

Jo z\z + X — 1 + ie) 



denote integrals over the pion distribution amplitude. Finally, q-- (xp) and Aq-- (xp) are the 
unpolarized and polarized valence and sea quark distribution functions evaluated at the 
point Xp ~ Q'^/{sx). The results for A, /i, and z/ within the pion bound state model are 
given elsewhere 0. 



III. NUMERICAL RESULTS 

The numerical results for the angular coefficient functions and u, eq. (|2.12|) , depend 
on X, p, the pion distribution amplitude and on the ratio of the polarized and unpolarized 
quark distribution functions of the proton. Whereas the unpolarized quark distributions 
are known quite well, the polarized ones still contain large uncertainties. For our results we 
choose the parametrizations given in at = 4 GeV^. 

In fig. 2 we present our results for several distribution amplitudes (f that are quite different 
in shape. Fig. 2a shows our choices for (p: The solid line represents the two-humped function 



p5| which gave a good fit to the data on A, fi and z/ ||^. The effective evolution parameter 
is set to ~ 4 GeV^. Since the predictions are also very sensitive to the behavior of the 



distribution amplitude in the endpoint region, we choose two different parametrizations ||T6 
for the convex distribution amplitude <f{z) = z"'{l — z)°'/B{a + l,a + 1). The dashed line 
in fig. 2a represents the asymptotical amplitude, i.e., a = 1, and the dotted line shows a 
narrow distribution amplitude with a = 10. 

In fig. 2b the moment / sin 26' sin 0(io"(s£ = +1) oc Np, is plotted in arbitrary units versus 
xl for the values s = 20^ GeV^, Q = 3 GeV and Qt = 0.9 GeV. It can be seen that 
different shapes of the distribution amplitude could be distinguished by a measurement of 
this quantity, which is proportional to the pion distribution amplitude evaluated at x ^ xl 
(cf. the numerator of ( p.l2D ). It is also demonstrated that narrow distribution amplitudes 
give a moment which vanishes in the large xl region. 

In figs. 2c and 2d the angular coefficient fi{se = +1) is plotted versus xl for two different 
choices of Qt and the same values for s and Q as in fig. 2b. In fig. 2c Qt was set again to 
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0.9 GeV, which gives the moderate value p = 0.3. It is demonstrated that the two-humped 
form for induces a minimum at xl ~ 0.6 which, however, vanishes for the smaller value 
p = 0.06 used in fig. 2d. For very narrow distribution amplitudes p is strongly suppressed. 
Since z/ = 2pp (cf. eq. ( p.l2| )) in the bound state model, we do not show separate plots for 



this quantity. 

Note that the effects of a variation of the pion distribution amplitude in the spin- 
dependent and spin-averaged cross sections - i.e. in the numerator and denominator of 
( p.l2| ) - partially compensate each other. As a consequence, the result for the two-humped 
distribution amplitude appears in between the results for the two convex distributions. The 
spin-dependent part of the cross section is a more sensitive "partonometer" , than the dimen- 
sionless coefficients p, v. Although the former quantity suffers from larger uncertainties due 
to higher order corrections, , one may expect only a small uncertainty for the determination 
of the shape of pion distribution amplitude, if the xl dependence of the corrections is weak. 

Although the naive parton model is at variance with the unpolarized angular distribu- 
tion we would like to present for completeness and comparison also the predictions of 
this model for p and v. We took the analytic results from and again used the parton 
distribution functions of [Q. In the range 0.5 < < 1 the parton model yields coefficients 
p{si = +1) and u{se = +1) which are positive and to a very good approximation indepen- 
dent of Xl. For the same values used in fig2b and 2c, i.e. s = 20^ GeV^,Q = 3 GeV and 
p = 0.3 we get a value of p{si = +1) ~ 0.12as ~ 0.036. This means that bound state effects 
are roughly of the same magnitude but opposite in sign compared to the parton model at 
moderate values of p. For the smaller value p = 0.06 we find p{se = +1) ~ 0.034as ~ 0.01, 
again independent of and positive. At such a small value of p the bound state model 
predicts that p and u are concentrated at large xl for not too narrow distribution ampli- 
tudes. In the kinematical range considered, the relation between u and p calculated within 
the parton model takes a similar form as in the bound state model, namely, u ~ 3pp. 

In conclusion we have shown that experiments on dimuon production from the scattering 
of pions on polarized targets will give detailed information both on the shape and on the 
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endpoint behavior of the pion distribution amphtude. 
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APPENDIX: RESULT FOR THE HADRONIC TENSOR 



Using the kinematic decomposition of the (symmetric part of the) hadronic tensor 

piipv pi^pi' ( piJ- pj^ A. pi^ pi") 



+ Qi seW,, 



(Al) 



where g'^" = g^^ — Q^Q^ /Q^ and = g^^P^ , we obtain from the bound state model the 
following predictions for the structure functions: 

W^i = ^ {([l - 5(1 + p')] F + (1 - x)Re I{x)f + 7rV(5)'| / {xp) 
W2 = -nx^f {xp) 

{(F + (1 - x)Re I{x)) [(1 - x{l + p^f)F + (1 - x)Re /(x))] + 7rV(5)^} 
1 - x(l + p')] ' F{F + Re I{x))f {xp) 
-p- l-x{l + F{F + (1 - x)Re I(x))f (xp) 



nx1,x 



1 — X 



~2 



1 — X 
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~3 



1 — a; 

2 ~2 

nXpX 



1 — X 



l-x(l + p2) 7r<^(x)FA/ (xp) . 



(A2) 



Here, 



n 



with Ci? 



2Nr 



is a normalization factor, and 



(A3) 
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9 
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(A4) 



are given by the unpolarized and polarized valence and sea quark distribution functions. 
The variable x and Xp as well as the integrals F and Re / containing the pion distribution 
amplitude are defined below eq. ( |2.12| ). 
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FIGURES 

FIG. 1. Leading contributions to the amplitude for the Drell-Yan process at large xl- 

FIG. 2. In (a) we show different pion distribution amplitudes ^{z). Solid line: 
two-humped function with an effective evolution parameter of ~ 4 GeV^; dashed line: 
ip{z) = — z)"/B(a + l,a + 1) with a = 1; dotted line: same with a = 10. In (b) the mo- 
ment / sin 20 sin <^c/o"(s^ = 1) is plotted versus xl for the values s = 20^ GeV^ Q = 3 GeV and 
Qt = 0.9 GeV. In (c) and (d) the angular coefficient fi{s£ = 1) is plotted versus xl for the two 
different values Qt/Q = 0.3 and Qt/Q = 0.06, respectively, and for the same values for s and Q 
as in (b). In (b-d) we use the same choices for (p{z) as in (a). 
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